######################
#  Replication code for 'Mediating the Electoral Connection', forthcoming in the JOP
#  John Henderson and John Brooks
#  12/7/2015    
######################    

# tableX-ghk.R
#  :: runs iv analysis to produce results for table X using gomez, hansford and krause rain measures

# ghk rain measure robust
  
rm(list=ls())
setwd('~/Dropbox/rainReplication')

fes.type=3
non.missings=4   

ghk=T

source('prelimRobust.R') 

#detach(covs)
# ANALYSIS             

#main_iv1_fe=ivreg(vote~
	#as.factor(fe_id_num)+
#	d_inc+#r_inc+
#	dist_prev + 
	#midterm + 
#	pres_party + 
    #black + #construction + educ + 
	#minc + farmer + forborn + gvtwkr + manuf + pop + unempld + 
	#urban + retail + 
	#sos + 
	#gov + 
	#comp_cq +
#	dose + dose_prv + vote_prv,
#	~#as.factor(fe_id_num)+
#	d_inc+#r_inc+
#	dist_prev + 
	#midterm + 
#	pres_party + 
    #black + #construction + educ + 
	#minc + farmer + forborn + gvtwkr + manuf + pop + unempld + 
	#urban + retail + 
	#sos + 
	#gov + 
	#comp_cq +
#	dose_prv + vote_prv + rain_ghk+rain_ghk_prev,
#	subset=full,data=covs)  
	  
main_iv2_fe=ivreg(vote~as.factor(fe_id_num)+d_inc+
	dist_prev + midterm + pres_party + 
	black + construction + educ + 
	minc + farmer + forborn + gvtwkr + manuf + pop + unempld + 
	urban + retail + sos + gov + comp_cq + redistricted + 
	dose + dose_prv + vote_prv,
	~as.factor(fe_id_num)+d_inc+
	dist_prev + midterm + pres_party + 
	black + construction + educ + 
	minc + farmer + forborn + gvtwkr + manuf + pop + unempld + 
	urban + retail + sos + gov + comp_cq + redistricted + 
	dose_prv + vote_prv + rain_elec00+rain_elec00_prev,
	subset=full,data=covs) 

#main_iv1_fe_sum=summary(main_iv1_fe)
main_iv2_fe_sum=summary(main_iv2_fe)
                                  
#main_iv1_fe_sumcl=coeftest.cluster(covs[full,],main_iv1_fe,cluster1='as.factor(fe_id_num)')
main_iv2_fe_sumcl=coeftest.cluster(covs[full,],main_iv2_fe,cluster1='as.factor(fe_id_num)')
      
save(
#main_iv1_fe,
main_iv2_fe,
#main_iv1_fe_sum,
main_iv2_fe_sum,
#main_iv1_fe_sumcl,
main_iv2_fe_sumcl,
file=paste('robustness/tableX_ghk-',non.missings,'_',fes.type,'.Rdata',sep=''))    

rm(
#main_iv1_fe,
main_iv2_fe,
#main_iv1_fe_sum,
main_iv2_fe_sum,
#main_iv1_fe_sumcl,
main_iv2_fe_sumcl
)
            
#detach(covs) 
#	}
#}    

# END      